Prediction of quantum stripe ordering in optical lattices 
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We predict the robust existence of a novel quantum orbital stripe order in the p-band Bose- 
Hubbard model of two-dimensional triangular optical lattices with cold bosonic atoms. An orbital 
angular momentum moment is formed on each site exhibiting a stripe order both in the superfluid 
and Mott-insulating phases. The stripe order spontaneously breaks time-reversal, lattice translation 
and rotation symmetries. In addition, it induces staggered plaquette bond currents in the superfluid 
phase. Possible signatures of this stripe order in the time of flight experiment are discussed. 
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Cold atomic systems with multiple components, such 
as large spin systems, exhibit much richer phase diagrams 
and properties than the usual spinless bosons and spin-i 
fermions. For example, various spinor condensations and 
spin dynamics have been investigated 1, 2\. Similarly, 
large spin fermions also exhibit novel features, includin 
hidden symmetries fij , quintet Cooper pairingstates 
and multiple-particle clustering instabilities [5j. Moti- 
vated by such considerations, but for orbital degeneracy, 
we investigate the p-band Bose- Hubbard (BH) model for 
cold atom optical lattices, finding a novel quantum or- 
bital stripe phase in triangular p-band optical lattices. 

In solid state physics, orbital dynamics plays impor- 
tant roles in transition metal oxides leading to inter- 
esting phenomena, such as orbital ordering and colossal 
magnetoresistence 0. In optical lattices, pioneering ex- 
periments on orbital physics have been recently carried 
out by Browaeys et. al [g by accelerating the lattice of 
bosons, and by Kohl et. al by using fermionic Fes- 
hbach resonance. These experiments demonstrate the 
population of higher orbital bands, motivating our theo- 
retical interest in possible orbital ordering in optical lat- 
tices. Compared to transition metal oxides where Jahn- 
Teller distortions often remove the orbital degeneracy, 
optical lattices have the advantage of the lattice rigidity, 
and thus the orbital degeneracy is robust. The p-band 
bosons in the sq uare or cubic lattices, have received much 
attention |^ 0, 0, 0| . For exaniple , R ef. focuses 
on the sub-extensive Z2 symmetry and the re- 

sulting nematic superfluid order by considering only the 
cr-type bonding in the band structure. By further keep- 
ing the TT-bonding term, the ground state is shown to 
break time reversal (TR) symmetry spontaneously, form- 
ing an antiferromagnetic order of orbital angular momen- 
tum (0AM) moments ITll.lT^. 

The p-band bosons in a frustrated optical lattice have, 
however, never been studied before. The experimental re- 
alization of the 2D triangular lattice has been discussed 



in the literature In this paper, we find a novel 

quantum stripe ordering of the p-band bosons in such 
a lattice. The onsite Hubbard interaction gives rise to a 
Hund's rule-like coupling in the 0AM channel, resulting 
in the formation of an Ising 0AM moment on each site. 
Due to the geometric frustration, the ground state ex- 
hibits a stripe order of the 0AM moments which sponta- 
neously breaks TR, lattice rotation and translation sym- 
metries. This stripe order bears some superficial similar- 
ity to its solid-state counterpart observed in strongly cor- 
related electronic systems, such as manganites JJ,], high 
Tc cuprates 0|, and high-Landau level quantum Hall 
systems but is different qualitatively since, unlike 
the solid-state examples, the stripe order in the p-band 
bosonic triangular optical lattices is fully quantum in na- 
ture and does not just arise from the long-range Coulomb 
interaction. 

We begin with the construction of the p-band BH 
model in a 2D (x-y) triangular lattice. The optical poten- 
tial on each site is approximated by a 3D anisotropic har- 
monic potential with frequencies uj^ uj^ — Uy = uJxy 
Thus we can neglect the px-band, and only consider a 
two-band model of p^; and py. We define three unit vec- 
tors ei = ix, 62,3 = "563; ± ^By, and the two primitive 
lattice vectors can be taken as aei.2 (a the lattice con- 
stant). The projection of the px,y orbitals along the 61,2,3 

directions are pi = Px,P2.3 — —-^Px ± -^Py Due to the 
anisotropic nature of the p-orbitals, the hopping terms 
are dominated by the "head to tail" type cr-bonding (the 
TT-bonding t± term can be neglected here) as 



Ho = U 



E { 

f,i=l,2,3 



pl,pPi,r+e 



h.c. 



(1) 



where is positive due to the odd parity of the p-orbitals. 
The on-site Hubbard repulsion Hint can be calculated 
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FIG. 1: The condensate configuration in the real space de- 
scribed by Eq.lPJ with TR, rotation, and U(l) symmetry 
breaking. The unit cell contains 4 sites as marked with thick 
lines. The OAM moments (dashed arrowed circles) form the 
stripe order, and induce bond currents exhibiting staggered 
plaquette moments (solid arrowed circles) . a in Eq. |3 de- 
pends on the interaction strength with q = ■|(-|) in the weak 
(strong) coupling limit respectively. Currents vanish in the 
weak coupling limit, and exist on tilted bonds at finite inter- 
action strength with directions specified by arrows. 

from the contact interaction of coupling constant g: 

r 

with the particle number and = ~i{p\py — p\px) 
the z-component OAM; U = /[A{2'iTfl'^ IJyl^] with 
li — ^Jfij (raixii) (i — x,y,z). The important feature of 
Hint is its ferro-orbital nature which is analogous to 
the Hund's rule for atomic electrons. This implies that 
bosons on each site prefer to go into the the axial states 
of px ± ipy This is because the axial states are spatially 
more extended than the polar states Px,y, and thus are 
energetically more favorable for (/ > 0. 

We first consider the weak couphng hmit, U/t — > 0. 
The Brillouin zone takes the shape of a regular hexagon 
with the edge length Air /(3a). The energy spectrum 

of Ho is E{k) = T ^/|-35fc}, where fj: = 

SLi cos(fc • ii) and 5^ = J23>i>j>i cos(fc • e^) cos(fc • Bj). 
The spectrum contains three degenerate minima located 
at = (0,^),if2,3 = (±f:;|z): The factor e'^-*' 
takes the value of ±1 uniformly in each horizontal row 
but alternating in adjacent rows. If the above pattern is 
rotated at angles of then we arrive at the patterns 
of e*^^'^ ''. Each eigenvector is a 2-component superpo- 
sition vector of px and py orbitals. The eigenvectors at 

energy minima are ipKi — s'^^^^lPy)- '/'-R'2,3 <^an be ob- 
tained by rotating ^pl at angles of ± ^ respectively. 

The ground state condensate wavefunction can be 
constructed as follows. Any linear superposition of the 
three band minima \l/c(r) = CiipKi + 02^/^X2 + c^i'K^ 
with the constraint |cip -I- |c2p + |c3p = 1 equally min- 



imizes the kinetic energy H^. However, an infinitesi- 
mal U jt removes the band degeneracy to further op- 
timize the interaction energy Hint- After straightfor- 
ward algebra, we find the optimal configurations occur at 
ci = 0, C2 = "^jCa = and its symmetrically equiv- 
alent partners. Thus the mean field condensate can be 

expressed as +'^'^k^] |0> with |0) the 

vacuum state and iVo the particle number in the conden- 
sate. This state breaks the t/(l) gauge symmetry, as well 
as TR and lattice rotation symmetries, thus the ground 
state manifold is U iV) ® Zi® Zj,. This state also breaks 
lattice translation symmetry, which is, however, equiva- 
lent to suitable combinations of C/(l) and lattice rotation 
operations. 

For better insight, we transform the above momentum 
space condensate to the real space. The orbital configu- 
ration on each site reads 

e*"^'"'(cosQ!|px) + icTf^sinalpj^)) (3) 

with a = ^ as J7/t ^ 0. The general configuration of a is 
depicted in Fig. Qfor later convenience. The J7(l) phase 

is specified at the right lobe of the p-orbital. The Ising 
variable ct? = ±1 denotes the direction of the OAM, and 
is represented by the anti-clockwise ( clockwise) arrow on 
each site. Each site exhibits a nonzero OAF moment and 
breaks TR symmetry. At \J jt — > 0, Px.y are not equally 
populated, and the moment per particle is This 
does not fully optimize Hint which requires iz,f — i^- 
However, it fully optimizes H^ which dominates over Hint 
in the weak coupling limit. We check that the phase dif- 
ference is zero along each bond, and thus no inter-site 
bond current exists. Interestingly, as depicted in Fig. ^ 
OAM moments form a stripe order along each horizon- 
tal row. The driving force for this stripe formation in 
the SF regime is the kinetic energy, i.e., the phase coher- 
ence between bosons in each site. By contrast, the stripe 
formation in high Tc cuprates is driven by the compe- 
tition between long range repulsion and the short range 
attraction in the interaction terms [l^. 

Next we discuss the ordering at large large values of 
V jt. We first minimize Hint with n particles per site. 
For simplicity, we consider the large n case, then Hund's 
rule coupling favors the onsite state | (cos \p\ + 

}n 
|0). This corresponds to the case a = ^ 

in Fig. [2 Because of the anisotropic orientation of the 
p-orbitals, the phase difference between two sites along 
each bond not only depends on the ?7(1) and the Ising 
variables, but also on the direction of the bond as in 
the p + ip Josephson junction arrays This effect 

can be captured by a U{1) gauge field. The effective 
Hamiltonian then reads iJoff = 12{ri r^) ^'^'^ K'^ri — 

- ^ri,r2 (cTfi , cTf^a)} + 3^ Si=-"-|) wheie the gauge field 
^n,r2 = ^yr^f^Ti.f^ — crf2^f2,ri; ^ri,f2 is the angle be- 
tween the bond from r\ to r2 and the x-axis, and thus 
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FIG. 2: A) Phase diagram based on the GMF theory in 
the 2x2 unit cell (see Fig. 0. Large scale GMF calcu- 
lations in a 30 X 30 lattice are performed to confirm the 
stripe ordered superfluid (SF) phase at points f, 2 and 3 
with {t/U,fi/U) = (0.02, 2), (0.03, f. 5) and (0.038,2.2), re- 
spectively. B) The flux $ around a rhombic plaquette v.s. 
t/(nU). It decays from i in the strong coupling limit to in 
the non-interaction limit. The solid line is the GMF result at 
n = 3, while the dashed line is based on the energy function 
Eq. Oof the trial condensate. 

^r2,fi = Ori,r2 + ^- The external gauge flux in the pla- 
quette i with three vertices ri.2,3 can be calculated as 

<I>i = -!- Ar^r' = ^(o-ri + CTfs + ^^^3 ) mod 1. (4) 

27r f-^ 6 

Following the analysis in Ref. [isl |20| , the ground state 
configuration for Ising variables requires the flux $i (vor- 
ticity) in each plaquette to be as small as possible, which 
are just ±i corresponding to Ising variables of two ±l's 
and one =Fl. 

The stripe order persists in the strong coupling SF 
regime due to the interaction among vortices. Consider 
a plaquette with vorticity -|-^, thus its three vertices are 
with two -l-l's and one —1. The neighbouring plaquette 
sharing the edge with two -|-l's must have the same vor- 
ticity, and merges with the former one to form a rhombic 
plaquette with the total vorticity +i. Thus the ground 
state should exhibit a staggered pattern of rhombic pla- 
quettes with vorticity of ±i. The stripe pattern of 0AM 
moments is the only possibility to satisfy this require- 
ment. The stripe phase obtained here is quite general: 
for example, it will also appear in the triangular-lattice 



p + ip Josephson-junction array if tunneling is dom- 
inated by the momentum-reversing process [20j, rather 
than by the momentum-conserving process that gives a 
uniform state. This stripe ordering possibility in the tri- 
angular Josephson junction arrays has not been ear- 
lier appreciated. 

In addition, the stripe order results in the staggered 
bond currents. We further optimize the U{1) phase vari- 
ables, and flnd that their pattern is the same as that in 
the weak coupling limit. Taking into account the equal 
weight of Px,y in each site, we find no phase mismatch 
along each horizontal bond, but a phase mismatch of 
A9 = ^ along each tilted bond. As a result, on each bond 
around the rhombic plaquette, the Josephson current is 
j = -^sinA^ where ng/n is the condensate fraction, 
and the current direction is specified by arrows in Fig. ^ 
The total phase winding around each rhombic plaquette 
is 4A6' — Itt, and thus agrees with the vorticity of ^. We 
emphasize that the p-band square lattice case does not 
have this interesting physics 1 1] . 

Since the stripe order exists in both strong and weak 
coupling limits, it should also exist at intermediate cou- 
pling strength. We have confirmed this conjecture using 
the Gutzwiller mean field (GMF) theory in a 30 x 30 lat- 
tice for three systems (marked as points 1,2, and 3 in Fig. 
Ejl. We find that the stripe ordered ground state is stable 
against small perturbations in all three cases. We further 
apply the GMF theory to the 2x2 unit ceU (Fig. ^1, and 
obtain the phase diagram of the stripe ordered SF and MI 
phases (Fig. [2K) . To understand the GMF numerical re- 
sults, we write the trial condensate with the p-orbital con- 
figuration on each site as e'^'^^ (cos a\px) + iapsmalpy)). 
It turns out that the pattern for the U{1) phase does 
not depend on a, and remains the same for all the cou- 
pling strength. The phase mismatch A9 on the tilted 
bonds reads A9 — 2^ — tt/2 with tan7 — v^tana, and 
the corresponding Josephson current is j = not sin A9. 
The value of a is determined by the minimization of the 
energy per particle of the trial condensate as 

£{a) ^ -t[\ + 2sin(2a + tt)] - -7- sin^ 2a -I- — . (5) 

D D o 

In the strong (weak) coupling limit, the energy minimum 
is located at a = ^ (|-), and thus the fiux in each rhom- 
bic plaquette $ = 4A6'/(27r) — (^), which agrees with 
the previous analyses. For the intermediate interaction, 
we present both results of $ at n = 3 based on the GMF 
theory and Eq.© in Fig.|33. They agree with each other 
very well, and confirm the validity of the trial condensate. 
Moreover, in the momentum space, the trial condensate 
for a general a can be expressed as -7== | (-0^_^ -|- 

|0), where ^'k,,,{^ = e'^^-^-i^^ ^^th 

\4>2,3,{ot)) = — cos.a\px) Tsinajpy) respectively. 

The QAM is of the Ising type, thus the stripe order- 
ing is robust against small perturbations such as a small 




(units of 2-nya) 

FIG. 3; Predicted TOF image for the stripe-ordered super- 
fluid phase with the condensation wavevectors K2,3- Note 
that the locations of the highest peaks depend on the size lx,y 
of the p-orbital Wannier function. Parameters are: a = 7r/6; 
lx,y/a — 0.1; the (5-function is replaced by a Lorentzian line 
for display. 

value of the 7r-type bonding t± <C . Further, we also 
check that the phase pattern in Fig. ^remains unchanged 
from minimizing the ground state energy. In particular, 
in the weak-coupling limit, the band minima and corre- 
sponding eigenvectors do not change at all in the pres- 
ence of a small t±, which renders the above conclusion 
obvious. 

The formation of the on-site orbital moment does not 
depend on the inter-site phase coherence, and thus the 
Ising variables can be ordered even in the MI state. We 
performed a ring exchange analysis showing the existence 
of the stripe-ordering at n > 2 provided the screening 
length of the interaction among vortices is larger com- 
pared to the size of the four-site plaquette. Due to loss 
of the inter-site phase coherence, bond currents disappear 
in the MI phase. 

Our predicted stripe phase should manifest itself in the 
time of flight (TOF) signal as depicted in Fig. O In the 
SF state, we assume the stripe ordering wavevector Ki, 
and the corresponding condensation wavevectors at K2,3. 
As a result, the TOF density peak position after a fight 
time of t is shifted from the reciprocal lattice vectors G 
as follows 

{n{r})t oc Y.[\Mc^,k)\^S^k-K2-G) 

G 

+ |(/>3(a,fc)|252(fc_ie3_G)}, (6) 

where k — mr/{ht); ^2,3 (a, fc) is the Fourier trans- 
form of the Wannier p-orbital wavefunction 1^2, 3(a)), fmd 
G — ^[to, (— m -I- 2n)/V3] with m, n integers. Thus 
Bragg peaks should occur at ^[m±i, -^(— m + 2n+ i)]. 
Due to the form factors of the p-wave Wannier orbit wave- 
function |(/)2.3(a, fc)P, the locations of the highest peaks 
is not located at the origin but around |fc| w l/lx,y Due 



to the breaking of lattice rotation symmetry, the pattern 
of Bragg peaks can be rotated at angles of In the 

MI phase, Bragg peaks disappear due to the loss of phase 
coherence. Instead, the stripe order appears in the noise 
correlations {n{r)n{r')) ^ which exhibit not only the usual 
peaks at G, but also peaks located at Ki + G. 

In summary, optical lattices provide a promising di- 
rection to study new phases of orbital physics. We focus 
on the p-band bosons in the frustrated triangular lattice 
which exhibits novel stripe orbital ordering of the on- 
site 0AM due to the geometric frustration effect to the 
phase coherence. In the SF phase, the staggered plaque- 
tte bond currents are also induced, reminiscent of the d- 
density wave proposal for the pseudogap phase in high Tc 
cuprates j2l| but with completely different microscopic 
origin. The stripe order persists in the MI phase even 
with the loss of superfluidity. The pattern of the stripe 
order can be observed in the TOF experiment. The or- 
bital physics of cold atoms opens up intriguing problems 
yet to be explored, for example, the p-orbital ordering in 
other frustrated lattices such as the Kagome. 
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